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Abstract 

A variant of the Vanderbilt ultrasoft pseudopotential scheme, where the 
normconservation is released for only one or a few angular channels, is pre- 
sented. Within this scheme some difficulties of the truly ultrasoft pseudopo- 
tentials are overcome without sacrificing the pseudopotential softness, i) 
Ghost states are easily avoided without including semicore shells, ii) The 
ultrasoft pseudo-charge-augmentation functions can be made more soft, iii) 
The number of nonlocal operators is reduced. The scheme will be most use- 
ful for transition metals, and the feasibility and accuracy of the scheme is 
demonstrated for the 4d transition metal rhodium. 



1 



Typeset using REVTgX 



With the development of the Vanderbilt Ultrasoft-pseudopotential (US) technique it 
has become possible to describe traditionally hard pseudopotential elements like transition 
metals and first row elements with a modest plane-wave cutoff. However, for some elements 
it is difficult to avoid the appearance of so-called ghost states without including semicore 
states in the pseudopotential construction 0. Furthermore, sometimes it is necessary to 
describe the US pseudo-augmentation charge with a higher cutoff than the pseudo-charge of 
the wave-function sum, and therefore so-called double-grid techniques have been developed 
0. Both these aspects significantly reduce the computational efficiency of the US scheme, 
and it is desirable to find a simple scheme to avoid these difficulties. 

In this report I will present a variant of the US technique, where only the Normconser- 
vation(NC) is released for some of the angular channels. The benefit from this construction 
is that the semi-local potential of one of the chemical active channels may be taken as the 
local potential. This reduces the number of projectors in the US scheme and by choosing the 
semi-local potential of the least bound reference state as the local potential, the appearance 
of ghost states is effectively avoided. Furthermore, since with this construction there are less 
pseudo-charge-augmentation functions, I have found the pseudo-charge smoothening proce- 
dure of Ref. || to be more efficient. In the following I will briefly review the main features 
of the US scheme, and present the modifications in a mixed US/NC scheme. Finally, I will 
build a US/NC pseudopotential for Rh and compare the transferability and computational 
efficiency with that of a NC pseudopotential and an US pseudopotential with semicore shells. 

In the US scheme the pseudo-wave-functions do not obey a NC condition and it is this 
feature which makes it possible to construct pseudopotentials with a modest plane-wave 
cutoff (< 30Ry) for traditionally hard pseudopotential elements. To make the pseudopoten- 
tial retain the first order scattering properties of the all-electron potential, as is the case of 
the NC pseudopotentials, the pseudo-wave-functions have to obey a generalized eigenvalue 
equation, where the missing norm enters in the overlap matrix. Furthermore, the density 
cannot simply be constructed from the pseudo-wave-functions, an augmentation charge has 
to be added in order to retain the correct electro-static potential in the interstitial region. 
Given a set of occupied pseudo- wave-functions the augmentation charge is given by 

r 9 (r) = E (*a|Xn>Qnm(r)<Xm|*«> , (1) 
a,n,m 

Q nm {r) = C(r)Var) - C(r)0 w (r) , (2) 

where ip n are the reference all-electron atomic wave-functions, <p n the corresponding reference 
atomic pseudo-wave-functions, and Xm are projectors upon the reference atomic pseudo- 
wave-functions. 

With this choice for the augmentation charge the total pseudo density does not only 
have the same norm as the all-electron density, as in the case of NC pseudopotentials, but 
is identical to the all-electron density in the case where the set of reference states inside 
each core form a complete basis for the all-electron wave-functions and the pseudo-wave- 
functions, i.e. ty a = X) n a nV'n an d $o = Hn^^n- To see this observe that since the pseudo 
and all-electron wave-functions and their radial derivatives coincide at the core radius they 
must have the same expansion in the reference states, thus a" = 6" which implies that 
T 1^ I 2 = T 1$ I 2 + o aug 

Therefore, in the US scheme not only the electro-static mono-pole of the core region is 
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correct, but also higher order poles are described rather accurately. Another feature of the 
US scheme is the possibility to include several reference states for each angular channel. 

In the NC scheme it is common practise to take one of the semi-local potentials as 
the local potential and thereby reduce the number of semi-local channels that have to be 
described. Furthermore, by choosing the semi-local potential of the least bound reference 
state one can effectively avoid the appearance of ghost states ||J| in a Kleinmann-Bylander 
U implementation. A similar choice for the local potential in the ultrasoft scheme would 
produce a very poor pseudopotential, since due to the missing norm of the corresponding 
pseudo-wave-function the local potential would only have the zero-order scattering proper- 
ties of the all-electron potential. Instead the common practise is to construct a semi-local 
potential for the first angular channel which is not chemical active and use that for the local 
potential, for instance, for the transition metals the semi-local potential of the /-channel is 
usually used for the local potential. However, for most transition metals this will produce 
a local potential with bound states far below the s- and p-reference eigenstates, and the 
appearance of ghost states is therefore unavoidable ||. A solution to this problem is to 
include s- and p-semicore shells in the pseudopotential construction, since these new refer- 
ence eigenstates usually are below the bound states of the local potential the ghost-states 
will disappear. The inclusion of semicore shells has the benefit that the accuracy of the 
pseudopotential gets comparable to all-electron calculations 0, however, it also increases 
the computer time and memory requirements of the calculation substantially. 

As an alternative I propose a mixed US/NC scheme, where the normconservation is only 
released for the hard channels, and the semi-local potential of one of the soft channels may 
thereby be taken as the local potential. For most transition metals the normconservation will 
only have to be released for the <i-channel, and the semi-local potential of the p-channel can 
therefore be taken as the local potential, thereby avoiding the appearance of ghost-states. 
The pseudo-charge-augmentation functions Q nm of Eq. (H) I now define as 



where £ = ip for the US channels and £ = <p f° r the NC channels, and only contributions 
from nonlocal channels are included, i.e. for most transition metals only the Qdd and Q s d 
elements will contribute to the augmentation charge. 

In the following I will construct a scalar-relativistic mixed US/NC pseudopotential for the 
4d transition metal Rh, show that the error introduced by the choice for the augmentation 
charge is minute, and compare the transferability and computational efficiency with a NC 
pseudopotential and a US pseudopotential with semicore states. 

To construct the NC pseudopotential (denoted NC(9)) I used the procedure suggested 
by Troullier-Martins M, with core radii 2.53, 2.53 and 1.39 (a.u.) for the s— , p- and d- 
channel, respectively. As suggested in Ref. 0, only the semi-local s— and d-potentials 
were constructed in the atomic ground-state configuration 4<i 8 5s 1 , while the semi-local p- 
potential was constructed in the 4d 7 5s°' 75 p ' 25 configuration. Figure 1 shows the pseudo 
wave-functions and their Fourier transforms. From the latter it can be seen that the s- 
and p-wave-functions are converged at ~ 20 — 25 Ry, while the d-wave-function must be 
described with a cutoff of ~ 50 — 60 Ry. From this pseudopotential the mixed US/NC 
pseudopotential (denoted US/NC (9)) was constructed by releasing the norm of the d- wave- 
function, increasing its core radius to 1.6 a.u, and then using the smoothening procedure 




(3) 
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of Ref. |10|j7[] to generate d-ultra-soft pseudo-wave-functions for two reference states, where 
one has the atomic eigenvalue (-0.4518 Ry) and the other the energy (0.2 Ry). With this 
construction the d-pseudo-wave-function can be described with a cutoff of ~ 25 — 30 Ry, as 
seen from its Fourier transform shown in Fig. lb(dotted line). Finally, I have constructed an 
US pseudopotential including 4s and Ap semi core shells (denoted US (17)), using 2 reference 
states for each channel, and a core-radius of 2.0, 2.0 and 1.6 for the s-, p-, and <i-channel, 
respectively. The US (17) pseudopotential is also converged at w 25 — 30 Ry. 

To generate pseudo- charge- augment at ion functions I have used the procedure described 
in Ref. 0,[7|, where the augmentation functions are replaced by L-dependent counterparts 

L 

nrai 

Qnm(r) = E cn LM Y LM(r)Q L nm (r), (4) 

LM 

which are smoothened inside a core radius, rf n , subject to the condition that the Lth moment 
of the electron charge density is conserved. One difficulty with this method is that if the 
core radius rf n is extended beyond a certain radius the smoothening procedure starts to 
develop negative sections in the Qn m ( r ) terms, leading to negative pseudo-charge densities. 
I have found that this difficulty is related to the different node structure of the radial all- 
electron wave functions, and that a larger r® n core radius can be chosen when there is only 
one angular channel. Figure 2 shows the Q® s , Qp p and Q° dd pseudo-charge-augmentation 
functions of the US (17) pseudopotential and the Q dd pseudo-charge-augmentation function 
of the US/NC(9) pseudopotential. For the US(17) pseudopotential the core radius could 
not be extended beyond rf n = 0.6 a.u., while the US/NC(9) core radius was r 4 ° n = 0.9 a.u. 
giving softer pseudo-charge-augmentation functions. 

Table 1 shows the transferability of the US/NC(9), NC(9) and US(17) pseudopotentials 
for atomic, bulk and surface properties of Rh. First notice the very high quality of the 
US (17) pseudopotential, having an accuracy essentially identical to all-electron calculations. 
Notice also the discrepancy (due to the local-density approximation) between the all-electron 
calculation and the experimental results for the bulk properties of Rh, which warns that the 
quality of a pseudopotential should never be judged by comparison with experimental data. 
The atomic calculations show that the main difference between the NC(9) and the US/NC(9) 
pseudopotentials, is the latter's enhanced description of the <i-electrons, which is due to the 
use of two reference states. 

Table lc shows the CPU time used to obtain the self-consistent charge density of the 
non-relaxed surface for the three different pseudopotential schemes. The comparison reveals 
that the US/NC(9) pseudopotential is 3.8 times faster that the US(17) pseudopotential and 
3.4 times faster than the NC(9) pseudopotential. For large systems the orthogonalization 
step is the most time-consuming part of a plane- wave program, and this scales as 0(NM 2 ), 
where N is the number of plane-waves and M the number of occupied bands. In the present 
case this gives a time factor relative to the US/NC(9) pseudopotential of (17/9) 2 ~ 3.6 for 
the US(17) pseudopotential and 2a ps 2.8 for the NC(9) pseudopotential. Another important 
issue is the memory requirement which scales as O(NM), which results in a factor 17/9 ~ 1.9 
for the US(9) pseudopotential and a factor 2a « 2.8 for the NC(9) pseudopotential relative 
to the US/NC(9) pseudopotential. 

To test whatever the neglect of the Q pd pseudo-charge-augmentation functions in the 
US/NC scheme introduce any significant error, I have constructed a US(17)* pseudopoten- 
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tial with Q sp = Q p d = Q s d = 0. These terms are most important for asymmetric systems, 
however, the calculated Rh(110) surface properties are almost unaffected by this approxi- 
mation(see Table lc). Furthermore, I have found that different variations of the US/NC 
scheme as including two s-reference states, using the atomic ground state for the p-reference 
configuration or including the Q ss term, do not change the quality of the pseudopoten- 
tial significantly. The discrepancy between the US (17) and US/NC(9) pseudopotentials is 
therefore mainly due to the latter's neglect of semicore states. 

In conclusion I have presented a simple scheme for constructing a mixed normcon- 
served/ultrasoft pseudopotential using as a starting point a normconserved pseudopoten- 
tial. The resulting pseudopotential scheme is more accurate and computationally efficient 
than the initial normconserved pseudopotential. Compared to truly ultrasoft schemes, ghost 
states are avoided without including semicore shells and the ultrasoft pseudo-augmentation 
charge is more easily described with the same cutoff as the wave-function pseudo-charge. 
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FIGURES 

FIG. 1. a) The real-space all-electron wave-functions of Rh and the pseudo- wave- functions of 
the NC(9) and US/NC(9) pseudopotentials. b) The Fourier-transformed pseudo-wave-functions. 

FIG. 2. The Fourier-transformed pseudo-charge- augmentation functions of the US (17) and the 
US/NC(9) pseudopotential. 
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TABLES 



TABLE I. Comparison of the three rhodium pseudopotentials. US(17) is the ultrasoft pseu- 
dopotential including semicore shells, NC(9) is the normconserved pseudopotential and US/NC(9) 
is the mixed normconserved/ultrasoft scheme. In all the calculations, including the atomic 
calculations of part a), the pseudopotentials are used in a non-local separable form, a) The 
calculated atomic-eigenvalues in different reference configurations and their deviation from the 
all-electron values, b) The lattice constant (ao), bulk modulus (B), derivative of the bulk mod- 
ulus (B'), and cohesive energy of rhodium, c) The surface energy (Esurf^) ? workfunction(TL), 
and interlayer relaxation (Ai2,A23 ) of the Rh(llO) surface. The US(17)* pseudopoten- 
tial neglects pseudo-charge-augmentation functions which couple different angular channels, i.e. 
Qsp — Qpd — Qsd — o. 

a) Atomic properties of Rh 

pseudo configuration 4d [Ry] A(4<i) [mRy] 5s [Ry] A(5s) [mRy] bp [Ry] A(5p) [mRy] 



US(17) 


Ad 9 bs°bp° 


-0.28267 


0.0 


-0.26864 


0.0 


-0.03501 


0.0 


US/NC(9) 


Ad 9 bs°bp° 


-0.28155 


1.2 


-0.27227 


-3.6 


-0.03798 


-3.0 


NC(9) 


4d 9 bs°bp° 


-0.27479 


7.9 


-0.27022 


-1.6 


-0.03883 


-3.8 


US(17) 


4d 7 5s 1 5p° 


-1.24280 


0.1 


-0.91783 


0.0 


-0.53999 


0.1 


US/NC(9) 


M 7 bs 1 bp° 


-1.24384 


-0.9 


-0.91472 


3.1 


-0.53760 


2.5 


NC(9) 


4d 7 bs 1 bp° 


-1.25787 


-15.0 


-0.91985 


-2.0 


-0.53800 


2.1 



b) Bulk properties of Rh 



pseudo 


cutoff [Ry] 


a [A] 


B [Mbar] 


B' 


E co /i[Ry/atom] 


Exp. 




3.80 a 


2.76 a 




0.42° 


All El. 




3.74 6 


3.46 b 


7.3 fc 


0.569 6 


US(17) 


30 


3.75 


3.44 


7.8 


0.636 


US/NC(9) 


30 


3.81 


3.16 


5.6 


0.609 


NC(9) 


60 


3.86 


2.98 


5.5 


0.594 



c) Surface properties of Rh(110) 



pseudo 


Cutoff [Ry] 


CPU time [s] 


E sur / [eV/atom] 


W [eV] 


Ad 12 [%] 


Ad 23 [%] 


Exp. 








4.98 c 


-6.9 d 


1.9 d 


US(17) 


30 


5796 


1.73 


4.77 


-12.1 


4.9 


US/NC(9) 


30 


1514 


1.92 


4.99 


-10.2 


2.5 


NC(9) 


60 


5162 


1.87 


4.98 


-10.5 


2.6 


US(17)* 


30 


5796 


1.73 


4.78 


-12.1 


4.9 



"Reference [11] ^Reference [12 c Reference JT5| (polycrystaline) ^Reference [14 
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